Autocorrelation exponent of conserved spin systems in the scaling regime following a 

critical quench 
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We study the autocorrelation function of a conserved spin system following a quench at the critical 
temperature. Defining the correlation length L{t) ~ t^'^, we find that for times t' and t satisfying 
L{t') <C L{t) <^ L{t')'^ well inside the scaling regime, the spin autocorrelation function behaves like 

(s(t)s(t')) ~ 1(1')-^''-^+'^^ [L{t')/L{t)]^'''. For the 0{n) model in the n ^ cx) limit, we show that 
Ac = d + 2 and — z/2. We give a heuristic argument suggesting that this result is in fact valid for 
any dimension d and spin vector dimension n. We present numerical simulations for the conserved 
Ising model in d = 1 and d = 2, which are fully consistent with the present theory. 



The quench of a ferromagnetic spin system [l|, from 
high temperature (To > Tc) to low temperature (typi- 
cally T = or T = Tc) is characterized by the growth of 
a correlation length scale (or domain length scale when 
domains can be identified), L{t) ~ i^/^. In the non con- 
served case, z depends on the final temperature of the 
quench (z = 2 for T < Tc, while z is the dynamical 
critical exponent for T = Tc j^l)- If the order parame- 
ter s(x,t) (possibly a vector) is locally conserved, z = 3 
(scalar) or z = 4 (vector) for a quench below Tc [l|, 
while z = 4 — ?7 J3 for a quench at Tc- Another inter- 
esting and fundamental quantity is the spin autocorre- 
lation A(i,0) = (s(x,i)s(x,0)) ~ L{t)-^ II mi- For 
non conserved dynamics, whatever the teinperature of 
the quench, A is non trivial (except in d = 1 [l|) and only 
approximate theories are available for T = [ij, 1^ , while 
for T — Tc j3j , the e-expansion of Ac can be calculated. In 
the case of conserved dynamics, it is now well established 
that Ac = A = d for quenches at and below Tc |y, L3> 1^ - 
Hence, for fixed t' and t -^ -|-oo, A{t,t') ~ L{t)-'^. How- 
ever, for t' and t > t' both in the scaling regime (in a 
sense to be defined later), several authors have observed 
numerically ^ .10. ill| and experimentally [l^l a faster 
power law decay of the autocorrelation. More precisely, 
in the case of a quench of an Ising system at Tc (critical 
quench), the authors of 11] obtained numerically 
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in d = 1 (where formally rj ~ I and Tc ~ 0) and d = 2. 
They respectively found AJ, ~ 2.5 in d = 1 and Ac ~ 3.5 
in d = 2. They also suggested a general scaling relation 



A(t,t') ^ Lity^c 
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where C{x) goes to a non zero constant for x 
C{x) - x-(^c-'*) for a; -^ 0, and 
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As noticed in [ll| , this scaling implies the existence of a 
new relevant length scale L{t')'^, which is the crossover 
length between the two observed regimes. Its physical 
meaning has yet to be elucidated. 

In the present Letter, we address the problem of the 
actual analytical derivation of Ac in the case of the 0{n) 
model in the limit of infinite n. Within this model, the 
diffusive nature of this new length scale can be under- 
stood, and one finds A'^ = d + 2 and = 2. By general- 
izing the interpretation of this diffusive crossover length 
scale to any 0{n) spin system, we conjecture that the 
result A^ = d -I- 2 holds and that (j> = 2 - 77/2 = z/2. 

We first examine the exactly solvable 0{n) model in 
the limit n —>■ 00 and for dimensions d > 2. This model is 
known to be pathological for a quench at zero tempera- 
ture, displaying multiscaling [13, whereas normal scaling 
should be restored at finite n [l|, |lj|. However, after a 
quench at Tc, the structure factor obeys standard scaling 
even for n ^ 00 [3 ■ In the standard Cahn-Hilliard equa- 
tion describing the evolution of the magnetization field 
s(x, i), s^(x, t)/n can be replaced by its average in the 
limit n — > 00. Thus, any spin component satisfies 
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where A:§ is a constant, r7(x, t) is a conserved delta- 
correlated noise satisfying {r]{\i.,t)r]{'k' ,t')) = 2Tcfc^(5(k-|- 
k')(5(i — i'), and (s^) has to be computed self-consistently. 
Although the derivation of the structure factor has al- 
ready appeared in the literature 7], we briefly repeat it 
as it furnishes a useful basis for our flnal derivation. 
Eq. 0) can be solved in Fourier space, leading to 



s(k,t) 
where 



s{k,0)+ [ e«('='^)77(k,r)dr 
Jo 



e-^C^'*', (5) 



q{k, t) = kH - P / [kl - (s2(x, r))] dr. (6) 

Jo 



Assuming an uncorrelated initial condition such that 
(s(— k, 0)s(k, 0)) = Sq, we then find the structure factor 

S'(A;,t) = (s(-k,i)s(k,i)) 



S{k,t) 



2Trk^ 
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-2q(k,t) 
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We now express the self-consistent condition (s^(x, t)) 



/ S{k,t) ,^^w , where A is the inverse of a lattice cut- 

oo) ~ fc~^+'', where rj is the 
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off. Tc is such that S{k,t 

usual critical exponent controlling the decay of the static 

correlation function (77 = for n -^ 00). This leads to 
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Finally, if the above condition is satisfied, we find that 
q{k, t) obeys a scaling relation for large t 



q{k, t) = q{kt^''^) = kh ~ Cdk^t^/"^, 



(9) 



where Cd is a universal constant determined by a simple 
integral relation {cd = for d > 4) J7|, and q{u) = u^ — 
CdV?- We thus find L{t) = t^/^, with 2; = 4, in agreement 
with the general result z = 4 — 77 ^3 • We hence reproduce 
the general form of the structure factor 



S[k, t) = s^e-^^t*^*'''' + t'^^-'^'^/'Fikt^^'-). (10) 

For the 0{n — 00) model, we have z = 4, ry = 0, and 



F{u) = 2Tcu' 



Note the following asymptotics for F{u) 
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In the scaling limit, the first term of the right hand side 
(RHS) of Eq. (fTn|l is negligible compared to the second 
term. In real space, Eq. pU|l illustrates the fact that 
conventional (critical) scaling is obeyed 

(s(x,t)s(0,t)) = L{tr('''^+^^)f[x/L{t)]., (14) 

where / is simply the inverse Fourier transform of F . 

We now move to the calculation of the two-time cor- 
relation function, focusing on the case where both con- 
sidered times t' and t > t' are in the scaling regime, a 
notion which will be made more precise hereafter. Us- 
ing Eq. 10), and working along the line of the derivation 
of S{k,t), we find the following expression for C{k,t) — 
(5(-k,i')s(k,i)) 



C{k,t,t') 
Ai{k,t,t') 

A2[k,t,t') 



Ai{k,t,t') + A2{k,t,t'), 

2 -q[kL{t')]-q[kL{t)] 

L{t'feiikL{t')]-qikLit)]p^,^m,'^^^ (17) 
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For a fixed t' and t ^ 00, the contribution of A2 becomes 
negligible, as for large t and hence L{t), only the contri- 
bution of small wave vector k ~ L{t)~^ matters. Using 
the result of Eq. p2|l , we indeed find that that this term 
is of order k'^ ~ L{t)^'^, whereas the main contribution 
Ai in Eq. (|16|) is of order Sq which is a constant. Contrary 
to what occurs for S{k,t), it is now the term depending 
on the initial conditions via Sq which dominates. Hence 
in this limit of fixed t' and t —> 00, we find 

C(fc, t, t') « C(fc, t, 0) = sge-«['=-^(*)l = G[kL{t)], (18) 

and in real space 

(s(x, i)s(0, t')) = m-^gix/Lit)], (19) 

where g is the inverse Fourier transform of G. One re- 
covers, in the limit t ^ i' to be made more precise later, 
that the large time autocorrelation exponent is Ac = d, 
which is observed in all conserved models including ther- 
mal fluctuations |y,l3- In this limit, conventional scaling 
holds. However, we will now show that the contribution 
of Eq. H17|l which has not so far been considered domi- 
nates in a well defined time regime, and will prompt us 
to introduce another autocorrelation exponent A^. 

For general t' and t > t', we now proceed to cal- 
culate the autocorrelation for a spin on a given lattice 
site. Defining A{t,t') = (s(x, i')s(x, i)), we finally find 
A{t,t') = Ai{t,t') + A2{t,t'), where 



Ai{t,t') 



-q[kL{t')]-q[kL{t)] 
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After a change of variable and noting that the region of 
k :» L{t)~^ barely contributes to the integral, we find 



A,{t,t') ^ L{t)-''a,[L{t')/L{t)] 



(21) 



Thus, Ai{t,t') obeys conventional scaling for any t' and 
t > t' . We explicitly find 



ai(-u) 



-k\l+u^)+cak-'{l+u^) 



d'^k 
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where this integral is now over the entire space. ai{u) 
remains bounded and of order Sq for any value of m = 
L{t')/L{t) < 1. Keeping the notation u = L{t')/L{t) < 
1, the expression for A2{t,t') can be written in the 
rescaled form 



A2{t,t') = L{t')-'^''-^^u'^x 

L(t)A 

^.k\l-u^)+c,k^-il^u-)p^^^^ 



(23) 
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Let us analyze the different asymptotics for ^2(^1 ^')- 
First of all, for large t — t' [u — \), the integral is 
dominated by the region of large fc's. Using Eq. ||HJ), 
we find the expected result A2(t, t) « /cq, which is the 



equilibrium value of (s^). Note that if i — i' ^ 1, 
we obtain A2{t,t') = A2(t,t) - Kd(t - t') + ..., where 
Kj^ is a computable constant. We now assume that 
1 ^ L{t) — L{t') <^ L{t'), which ensure that u is very 
close to 1. In this regime, we find that 



A2{t,t') - JdL{t) 
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, (24) 
(25) 



where Jd and J'd can be written exactly as simple inte- 
grals. Finally, and this constitutes the central result of 
this Letter, we consider the regime 1 <C L{t') <C L{t). In 
this case, u 4C 1, and the integral of Eq. (|23|l is dominated 
by the region of k of order unity, so that the small ar- 
gument asymptotics can be taken for F{ku) in Eq. (|23|l . 
We find 



A2{t,t') ^ KdL{t')-'^^-^^ 



L[t') 
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(26) 
(27) 



Eq. (|26|l takes exactly the expected form of Eq. iQJ , with 



Xl. = d + 2. 



(28) 



Hence, we find that A2{t,t') prevails over Ai{t,t') for 
L{t') < L{t) <C Lo(t'), with Lo(i') ^ L{t')'^ and = 2. 
For L(t') < L(t) < Lo(t'), the autocorrelation A{t, t') k, 
A2{t,t') is then given by Eq. ^. Moreover, Eq. ^ 
shows that instead of Eq. Q , the correct scaling is rather 

A{t, t') = A{t, 0) + L{t')-'^'^-^+'^^D[L{t)/L{t')], (29) 

with D{l/u) ^ u^" for u ^ 1. Both scaling are equiva- 
lent only for w <C 1. We now present an heuristic argu- 
ment based on dimensional analysis which suggests that 
the result \'^ = d + 2 may be of general validity for con- 
served spin systems. Indeed, the occurrence of a new 
length scale bigger than L{t) could have been inferred 
from the small k behavior of S{k, t). In the n ^ oo limit 
and for fc -+ 0, Eq. ^^ leads to 



S{k,t) 



2Tck^L(t)'^ 



(30) 



A natural momentum scale fco(^) '^ io(^)^^ arises by 
matching the two terms of the RHS of Eq. H3U|) . which 
leads to = 2 and hence A^ = d -|- 2. 

In the general case, for short-range correlated initial 
conditions, we expect the following general form to hold 

S{k,t) ^ Fi[kL{t)] + Lit)^-'^F2[kLit)], (31) 

with Fi(0) = Sq being a non zero constant (equal to the 
variance of the initial total magnetization normalized by 
the volume), while the scaling contribution should vanish 
for fc = 0, implying -F2(0) = 0. Imposing F2{p) ~ p'', 7 



< 10-' - 




L(t)/L(t') 



FIG. 1: Illustrating the result of Eq. ^, we plot A{t,tk) - 
yl(i,0) = D[L(t)/L(tfc)], for L(t^) ^ L{0)r^ , with r = 1.75 
and k = 1,...,7 (40000 samples of length A^ = 5000). Al- 
though the initial slope is smaller (Ac ~ 2.5 [ij; dotted line), 
the asymptotic exponent is very close to A^ = 3 (dashed line 
fit). The bottom inset shows the original data for A{t,tk) 
and A{t,tk) — A{t,0) (dashed lines). The top inset shows 
L{t)[A{t,tk) - A{t,0)] as a function of L{t)/L{tk)'^ (with 
cj> = 3/2). Lines of slope A^ — 1 are shown for AJ, = 3 (dashed 
line) and A^ « 2.5 (dotted line). 



is necessarily an even integer. If 7 were not integer, the 
correlation function scaling function / defined in Eq. I|14|) 
would have a power law decay for large distance, which is 
unphysical as such correlations cannot develop in a finite 
time starting from short-range ones. 7 cannot be an odd 
integer as space isotropy guarantees that / should be an 
even function. Contrary to the case of a quench at T = 
0, for which convincing theoretical arguments for d > 2 
[l5j and experiments llQ] show that F2{p) ~ p'^, there 
is no reason to expect the same for critical quenches. 
Generically, we expect F2(j>) ~ p'^ as found for the d = 1 
conserved Ising model |a, LSI i and in the present Letter for 
the 0{n) model for n ^ oo. Finally, the small k behavior 
of the structure factor should be of the form 



S{k,t) 



Cak'Litf-^ + 



(32) 



where Co > is a constant. Assuming that the length 
scale obtained by matching both terms of the RHS of 
Eq. (|32|l is the same as the crossover length between the 
two observed regimes for the autocorrelation, and using 
the general result of Eq. Q , we obtain 



2 - 7?/2 = 1 



V 



X'-d' 



(33) 



which implies A^ = d + 2. This result also extends to 
d ~ 1 after formally taking rj — 1, leading to A^ = 3 
and (p ~ 3/2. Note that the crossover scale can also be 
written Lo{t) ~ i*/^ ~ i^/^, which is the diffusion scale 
associated to thermal noise. At least in d = 1, this scale 
can be related to the equilibrium diffusion of tagged spins 
observed in |l3] ■ We now present simulations of the Ising 
model Kawasaki dynamics in d = 1 and d = 2 after a 




L(t)/L(t') 



FIG. 2: In the bottom inset, we plot A{t,tk) and A{t,tk) — 
A{t,0) (dashed hnes) as a function of L{t), for L{tk) ~ r'' , 
with r = 1.5 and k — 1, ..., 5 (16 samples of size N = 500x500, 
L(0) = 1/^2). The main plot shows L{tk)'^[A{t,tk) - 
A{t,0)] = D[L{t)/L{tk)]. Although the initial slope is con- 
sistent with A^ ~ 3.5 11] (dotted line), the effective exponent 
certainly increases and the asymptotic slope is more com- 
patible with Ac = 4 (dashed line fit). The top inset shows 
L{t)'^[A{t,tk) - A{t,0)] as a function of L{t)/ L{tk)''' (with 
= 15/8). Although not as clean as in d = 1, the scaling 
plot is better described by the line corresponding to A^ = 4 
(dashed line) rather than A^ ~ 3.5 (dotted line). 



quench at Tc- In the d — 1 case, we use the accelerated 
algorithm introduced in [7| , which is faster than that used 
in lljj (but does not permit to compute simply the re- 
sponse function as was needed in [ll|). By fitting A{t, t') 
in the scaling regime, the authors of [Il| found A^ ~ 2.5 
lower than our prediction A^ = 3. However, for the mod- 
erately large numerically accessible times, the contribu- 
tion of Ai{t,t') w Ait^Q) is significant. When plotting 
A{t, t') — A{t, 0) as a function of L{t), one actually finds 
AJ, w 3 instead of A'^ « 2.5. Result of simulations for 
the d = 2 Ising model evolving with Kawasaki dynamics 
at Tc are shown on Fig. 2. Considering the very slow 
growth of L(i) ~ i'*/^^, it is difficult to obtain data span- 
ning more than one decade in L{t). Hence, the regime 
of interest 1 <^ L{t') <^ L{t) cannot be reached and the 
separation of scales properly achieved. Still, subtracting 
A(t, 0) from A{t,t') leads to A^ sa 4, significantly greater 
than the value A^ « 3.5 found in [ll|. 

In conclusion, in view of the exact result for the 
0{n = oo) model, a general argument for any n and 
d, and convincing simulations in d = 1 (and consistent in 
d = 2), we have strongly suggested that A^ = d -f 2 and 
(f) — z/2 generally holds. We also find that the scaling 
form of Eq. (|29|l is more appropriate than Eq. |(2Jl. The 



compelling generalization of our heuristic argument to a 
quench at T < Tc (in d > 2, and admitting F2{p) ^ p*) 
leads to A{t,t') - [L{t')/L{t)]^' for L{t') < L{t) <C 
L(t')'^, with A' = d + 4 and (/) = 1 -f d/4. In d = 2, 
the prediction A' = 6 is significantly larger than the nu- 
merical result A' « 4 J9(. However, the fit in J9( was 
performed in the short scaling regime over less than a 
decade in L{t), and subtracting A{t,0) before perform- 
ing the fit could lead to a significantly higher value for 
A', as noted in the two examples treated in this Letter. 
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